Mode density of orbital angular momentum modes in a cylindrical cavity and in free 
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In this short note we derive an approximate expression for the mode density of modes carrying 
orbital angular momentum (OAM) in a cylindrical cavity which is large compared to the inverse 
wavevector in any direction. We argue that in the large cavity limit the modes asymptotically 
converges to the corresponding OAM modes in free space. We compare the result to Weyl’s law. It 
is found, as expected, that the mode density does not depend on whether or not the modes carry 
OAM. 

I. INTRODUCTION 

Recently, the interest in using the angular momentum degree of freedom in both classical and quantum optics has 
increased. The attraction for both communities is the anticipated possibility to utilize angular momentum as a new 


>-^degree of freedom. Claims have been made that using the orbital angular momentum degree of freedom will increase the 
' ^ communication capacity per unit frequency and unit volume in radio links [l|. We cite from the abstract of [l|: “This novel 
radio technique allows the implementation of, in principle, an infinite number of channels in a given, fixed bandwidth, even 
' without using polarization, multiport or dense coding techniques.” This claim has generated considerable excitement in 


g the physics, electromagnetics, and microwave engineering connnnnities m -d therefore it merits a critical examination, 

Q\ Weyl’s law states that any set of orthogonal modes have the same density of states provided the “cavity” or volume 
^ one confine the modes into is large enough. In fact, the requirement is that for any spatial co-ordinate q, the cavity 
(N dimension q in that direction times the wave vector in the same direction kq should fulfil qkq 1. “Free space” certainly 


in 


fulfills the condition as one may make a fictions cavity arbitrarily large and let it have an arbitrary shape. 


To exemplify Weyl’s law we derive the mode density for two sets of modes. In Sec. Ullwe derive the mode density for a 
. ^ cubic cavity. Although the final result is very well known from the literature we do this by splitting the density of states 
into a transverse and a longitudinal factor. Such a division simplifies the corresponding derivation in Sec. IIIIl where we 
^ derive the mode density of angular momentum modes in a cylindrical cavity and show that the mode density shows no 
increase. The latter modes are essentially equal the OAM Bessel modes in free space, so the claim of an communication 
capacity increase in free space is unfounded. 


II. DENSITY OF STATES IN A CUBIC CAVITY 

To quickly derive a frame of reference without consideration of angular momentum, we start by calculating the density 
of modes in an empty, ideal, metallic cavity formed like a cube with the side L. The modes in such a cavity are TE and 


2 


TM modes satisfying the dispersion relation 



( 1 ) 


where uj is the angular frequency, c is the speed of light, and m,n, and q are 0, 1, 2, ... but two of them cannot be zero 
at the same time. Hence, along all axes in cartesian fc-space the mode spacing is tt/L. 

Now let us divide the A:-vector into a transverse and a longitudinal direction. We will arbitrarily take kz as the 
longitudinal wavevector component and the component lying in the ccy-plane to be the transverse component k±. The 
number of TE-modes Nj_ having a transverse wavevector smaller or equal to k± is found by dividing the circular k-space 
area irk^ with the transverse fc-vector “unit cell” (tt/L)^, that is: 


N± 




( 2 ) 


Hence, 


dN± _ 2L^k± 
dk± IT 


( 3 ) 


The mode density in the longitudinal direction is dNz/dkz = Lj-K. To get the total number N (i.e., both TE and TM) 
of modes with a wavevector shorter or equal to k we integrate the mode density over all such wave vectors. We keep in 
mind that k'^ = k‘]_ + k^. 


N = 2- 
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where the pre-factor 2 in front of the first integral accounts for the inclusion of TM modes, and the subsequent factor 
1/8 is due to that only modes in the first octant of /c-space should be counted since m, n, and q are positive integers. 
Thus, we get the mode density 


dN 

dk 


( 5 ) 


where V = is the cavity volume. This density can trivially be recast as 

dN _ SttVv^ dN _ Vuj^ 
dv du) ’ 


where v denotes frequency, w denotes angular frequency, and c is the speed of light in vacuum. The only approximation 
in JS]) is due to the discrete nature of m, n, and q when the cavity side length L is on the order of a wavelength (so that 
Vk^ is on the order of unity). When 1, ([S]) is exact. 


III. ORBITAL ANGULAR MOMENTUM MODES IN A CAVITY WITH CIRCULAR CROSS-SECTION 

The modes inside a hollow, circular, metallic waveguide can carry angular momentum in the form of orbital angular 
momentum and spin. In many ways they are very similar to the modes of a Bessel beam (that also can carry 0AM) 
and when the cavity dimensions are large compared to cjv the solutions are essentially identical to Bessel beams in free 
space Q. Therefore, under this condition the mode density is essentially equal for the free-space modes and the cavity 
modes. However, we find that the cavity mode density is easier to derive which is why we look at this case. 
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In a circular waveguide with radius a the TE modes are given by the expressions (in cylindrical coordinates) 

_ iAk^ / Pi^p \ 

\ a J 

^ _ iAkf^l / Pl^p \ ^{ilip) ^(i[ujt-Pir„,z]) ^ 

^ImP V “ / 
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( 7 ) 
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where A is an amplitude (with SI unit Volt), po denotes the permeability of vacuum, J((a:) is the first derivative of the 
hth Bessel function of the first kind, gives the m:th solution to J/(x) = 0, kim = p\ml^ ^^.d k^ = oj/c. One can show 
that these eigenmodes carry the angular momentum Ih per photon, where I is an integer (but can take on both positive 
and negative values) and furthermore, when ^ 1 the angular momentum is primarily carried as orbital angular 

momentum [3 • It is important to remark that if one pair-wise superimposes such modes with equal amplitude and radial 
mode index but opposite chirality (Z and —Z), then these new modes carry no 0AM. The in-phase superposition will 
furthermore be orthogonal to the out-of-phase superposition so we will get an OAM-free, complete set of modes with 
pairwise correspondence to the 0AM modes (except for the Z = 0 mode that has no 0AM anyhow). 

If we make a cavity of such a waveguide of length L by sealing both ends with a perfectly conducting plate, then in 
the longitudinal direction the cavity eigensolutions will be standing waves with a longitudinal mode density L/tt. To 
count the number of modes in the radial direction we will make some approximations of the Bessel functions of the first 
kind, since the functions and their zeros are not easily expressible analytically. First we note that Ji{x) and J/(a;) are 
oscillating around zero with roughly the period tt for large values of x (this is a slight underestimation of the period). This 
means that the spacing between consecutive p'l^ (and also the corresponding zeros for Ji(x)) for a given Z must also be 
very nearly tt. This gives us the mode spacing between consecutive transverse wavevectors, e.g., — p'lm)!^ = t: ja 

which translates to a mode density of o/tt. 

Secondly, the first zero of Ji{x) and J/(cc) occurs roughly when x = 1. For example, the first zero of J 2 ooo(a^) occurs 
for X K, 2023,46 and the first zero of J 2 ooi(a^) occurs for X k, 2024,47. Thus, with an increment of 1/a in wavevector 
space, one will exceed the wavevector cut-off for another series of mode solutions with even higher angular momentum. 
Summing the number of modes with a transverse wavevector no larger than fc_L one gets (see Fig. [Ij 
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FIG. 1: The cylindrical cavity modes marked in a figure where the horizontal coordinate gives the transverse wavevector k± and 
the vertical coordinate gives the OAM mode index 1. 
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N o fej 2a^fci 

^ (7r/a)(l/a) tt ’ 

where the pre-factor 2 again takes into account that there is also an equally large set of TM solutions 
transverse mode density is 

dNj_ 4a^kj_ 

dk± TT 


The total number of modes are thus 


4a^L 

N = —^ / dkz 


4^a^Lk^ 


Taking the cavity volume V = na^L into account, the mode density approximated this way becomes 


dN 

dk 


4Vk^ 


4 Vk^ 


(9) 


Hence, the 


( 10 ) 


( 11 ) 


( 12 ) 


We see that except for the pre-factor 4/7r si 1.27 we arrive at the same result as from the previous derivation. The 
discrepancy between m and uni) is not surprising, because we have used approximations to derive the mode spacing for 
these angular momentum modes whereas the derivation of (jS]) was exact (under the condition stated in Sec. |TT|). That our 
estimate is an overestimate is expected, because the distance between successive zeros for J/(a;) is in fact always larger 
than TT, and as we have seen above, the first zero of the function occurs not when x = I, but for a somewhat larger value. 
Therefore we have underestimated the mode spacing leading to an overestimation of the mode density. This is reflected 
by the numerical pre-factor 4/7r that is an artefact due to our approximations. 


IV. CONCLUSIONS 

We have derived an approximate expression for the mode density for a set of complete, orthogonal modes in a cavity 
with a circular cross section. These modes carry angular momentum and the mode index I directly gives the number 
of angular momentum quanta per photon the modes carry. As explained above, this mode set therefore very closely 
approximates free-space Bessel beams when the cavity dimensions are much larger than cjv and thus should have 
the same mode density. We find that the mode density of these modes has the same scaling with frequency, angular 
frequency or wave vector as, for example, a plane wave set of modes would have. From this we conclude that the angular 
momentum carried by electro-magnetic waves does not increase the channel capacity compared to a mode set void of 
orbital angular momentum. (Such a mode set could easily be constructed from equally weighted sums and differences of 
angular momentum modes with indices I and —I [^). Using one set of eigenmodes one codes the information in another 
degree of freedom than with the other set, but the number of power-orthogonal modes (or “channels”) per unit frequency 
bandwidth remains the same. 
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